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We present calculations of the exchange interactions and Curie temperatures in Cr-based pnictides 
and chalcogenides of the form CrX with X=As, Sb, S, Se and Te, and the mixed alloys CrAssoXso 
with X=Sb, S, Se, and Te. The calculations are performed for Zinc Blende (ZB) structure for 12 
values of the lattice parameter between 5.44 and 6.62 A, appropriate for some typical II- VI and 
III-V semiconducting substrates. Electronic structure is calculated via the linear muffin-tin-orbitals 
(LMTO) method in the atomic sphere approximation (ASA), using empty spheres to optimize ASA- 
related errors. Whenever necessary, the results have been verified using the full-potential version of 
the method, FP-LMTO. The disorder effect in the As-sublattice for CrAs 5 oX 5 o (X=Sb, S, Se, Te) 
alloys is taken into account via the coherent potential approximation (CPA). Exchange interactions 
are calculated using the linear response method for the ferromagnetic (FM) reference states of the 
alloys, as well as the disordered local moments (DLM) states. These results are then used to estimate 
the Curie temperature from the low and high temperature side of the ferromagnetic/paramagnetic 
transition. Estimates of the Curie temperature are provided, based on the mean field and the more 
accurate random phase approximations. Dominant antiferromagnetic exchange interactions for some 
low values of the lattice parameter for the FM reference states in CrS, CrSe and CrTe prompted us 
to look for antiferromagnetic (AFM) configurations for these systems with energies lower than the 
corresponding FM and DLM values. Results for a limited number of such AFM calculations are 
discussed, identifying the AFM[111] state as a likely candidate for the ground state for these cases. 
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I. INTRODUCTION 

Half-metallic ferromagnets with high (room temper- 
ature and above) Curie temperatures T c are ideal for 
spintronics applications, and as such, much experimen- 
tal and theoretical^— effort has been devoted in recent 
years to the designing of and search for such materials. 
Among these, Cr-doped dilute magnetic semiconductors 
(DMS^ or Cr-based alloys and in particular CrAs and 
CrSb£~— in Zinc Blende(ZB) structure have attracted 
particular attention, not only because of the possibility 
of complete spin polarization of the carriers at the Fermi 
level, but also for their possible high T c . Akinaga et al~ 
were able to grow ZB thin films of CrAs on GaAs (001) 
substrates by molecular beam epitaxy, which showed fer- 
romagnetic behavior at temperatures in excess of 400 K 
and magnetic moments of 3/j.b per CrAs unit. Theoret- 
ical calculations by Akinaga et al& and several other 
theoretical calculations since ther£ ' 12 ' 13 ' 17 ~ — have veri- 
fied the half-metallic character of CrAs. The high value of 
T c has also been supported by some of these studies»iZr— 
Thin films of CrSb grown by solid-source molecular beam 
epitaxy on GaAs, (Al,Ga)Sb, and GaSb have been found 
to be of ZB structure and ferromagnetic with T c higher 
than 400 

Galanakis and Mavropoulos,— motivated by the suc- 
cessful fabrication of ZB CrAs, CrSb and MnAs^i have 
examined the possibility of half-metallic behavior in or- 



dered ZB compounds of transition metals V, Cr and Mn 
with the 'sp' elements N, P, As, Sb, S, Se and Te. Their 
theoretical study shows that the half-metallic ferromag- 
netic character of these compounds is preserved over a 
wide range of lattice parameters. They also found that 
the half-metallic character is maintained for the transi- 
tion element terminated (001) surfaces of these systems. 
Yamana et alJ£ have studied the effects of tetragonal 
distortion on ZB CrAs and CrSb and found the half- 
metallicity to survive large tetragonal distortions. Of 
course, the ground states of many of these compounds in 
the bulk are known to be different from the ZB structure, 
the most common structure being the hexagonal NiAs- 
type. Zhao and Zunger— have argued that ZB MnAs, 
CrAs, CrSb, and CrTe are epitaxially unstable against 
the NiAs structure, and ZB CrSe is epitaxially stable 
only for lattice constants higher than 6.2 A, remaining 
half-metallic at such volumes. They also find that even 
though the ground state of CrS is ZB, it is antiferro- 
magnetic at equilibrium lattice parameter, and thus not 
half-metallic. These results reveal the challenge exper- 
imentalists face in synthesizing these compounds in ZB 
structure. However, the possibility remains open that 
such difficulties will be overcome with progress in tech- 
niques of film-growth and materials preparation in gen- 
eral. Recently, Deng et al^ were successful in increasing 
the thickness of ZB-CrSb films to ~ 3 nm by molecular 
beam epitaxy using (In,Ga)As buffer layers, and Li et 
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ali^ were able to grow ~ 4 nm thick ZB-CrSb films on 
NaCl (100) substrates. 

In view of the above situation regarding the state of ex- 
perimental fabrication of these compounds and available 
theoretical results, it would be appropriate to study the 
variation of magnetic properties, particularly exchange 
interactions and the Curie temperature, of Cr-based pnic- 
tides and chalcogenides as a function of the lattice pa- 
rameter. Towards this goal, we have carried out such 
calculations for the compounds CrX (X=As, Sb, S, Se 
and Te) and the mixed alloys CrAssoXso with X=Sb, S, 
Se and Te. Essentially we study the effect of anion doping 
by choosing elements of similar atomic sizes (neighboring 
elements in the Periodic Table), one of which, namely Sb, 
is isoelectronic to As, while the others (S,Se,Te) bring 
one more valence electron to the system. The mixed 
pnictide-chalcogenide systems offer further opportunity 
to study the effects of anion doping. The alloying with 
other 3d transition metals (both magnetic, e.g. Fe or Mn, 
or non-magnetic , e.g. V) on cation sublattice would also 
change the carrier concentration and bring about strong 
d-disorder which can additionally modify the shape of 
the Fermi surface. This, however, is not the subject of 
the present paper. 

Almost all theoretical studies on these alloys so far 
address aspects of electronic structure and stability of 
these alloys only. Although a few theoretical estimates 
of exchange interactions and the Curie temperature for 
CrAs at equilibrium lattice parameter have appeared in 
the literature, a detailed study of the volume depen- 
dence of these quantities is missing. For the other alloys, 
CrSb, CrS, CrSe and CrTe, no theoretical results for the 
exchange interaction, Curie temperature and their vol- 
ume dependence exist at present. The mixed pnictide- 
chalcogenide systems offer the possibility of not only 
creating these alloys over a larger range of the lattice 
parameter, but also with a larger variation in the ex- 
change interactions. This is because at low values of the 
lattice parameter the dominant Cr-Cr exchange interac- 
tions in the chalcogenides can be antiferromagnetic, while 
for the pnictides they are ferromagnetic. The pnictide- 
chalcogenide alloying is important from the experimental 
viewpoint of stabilizing the ZB structure on a given sub- 
strate, via the matching of the lattice parameter of the 
film with that of the latter. Although the present study 
is confined to the ZB structure only, we hope that it will 
provide some guidance to the experimentalists in their 
search and growth of materials suitable for spintronic de- 
vices. 



II. ELECTRONIC STRUCTURE 

Electronic and magnetic properties of CrX (X=As, 
Sb, S, Se and Te) and CrAs 5 oX 5 o (X= Sb, S, Se and 
Te) were calculated for lattice parameters varying be- 
tween 5.45 and 6.6 A, appropriate for some typical II- 
VI and III-V semiconducting substrates. Calculations 



were performed using the TB-LMTO-CPA method 2 ^ 
and the exchange-correlation potential given by Vosko, 
Wilk and Nusair^ In our LMTO calculation we opti- 
mize the ASA (atomic sphere approximation) errors by 
including empty spheres in the unit cell. We use the fee 
unit cell, with Cr and X (As, Sb, S, Se and Te) atoms 
located at (0,0,0) and (0.25,0.25,0.25), respectively, and 
empty spheres at locations (0.5,0.5,0.5) and (-0.25,-0.25,- 
0.25). For several cases, we have checked the accuracy 
of the LMTO-ASA electronic structures against the full- 
potential LMTO results^ and found them to be satis- 
factory. For the mixed alloys CrAs5oX 50 (X= Sb, S, Se 
and Te), the As-sublattice of the ZB CrAs structure is 
assumed to be randomly occupied by equal concentra- 
tion of As and X atoms. The disorder in this sublattice 
is treated under the coherent potential approximation 
(CPA). 

Our spin-polarized calculations assume a collinear 
magnetic model. In the following we will present re- 
sults referred to as FM and DLM. The FM results fol- 
low from the usual spin-polarized calculations, where self- 
consistency of charge- and spin-density yields a nonzero 
magnetization per unit cell. Although we call this the 
FM result, our procedure does not guarantee that the 
true ground state of the system is ferromagnetic, with the 
magnetic moments of all the unit cells perfectly aligned. 
This is because we have not explored non-collinear mag- 
netic states, nor all antiferromagnetic (AFM) states at- 
tainable within the collinear model. Indeed, our results 
for the exchange interactions in some cases do suggest the 
ground states being of AFM or complex magnetic nature. 
For lack of a suitable label, we refer to all spin-polarized 
calculations giving a nonzero local moment as FM state 
calculations. Within the Stoner model, a nonmagnetic 
state above the Curie temperature T c would be character- 
ized by the vanishing of the local moments in magnitude. 
It is well-known and universally accepted that the neglect 
of the transversal spin fluctuations in the Stoner model 
leads to an unphysical picture of the nonmagnetic state 
and a gross overestimate of T c . An alternate description 
of the nonmagnetic state is provided by the disordered 
local moment (DLM) model, where the local moments 
remain nonzero in magnitude above T c , but disorder in 
magnitude as well as their direction above T c causes the 
global magnetic moment to vanish. Combining aspects of 
the Stoner model and an itinerant Heisenberg-like model, 
Heine and co-worker a 28 ' 29 have developed a suitable crite- 
rion for a DLM state to be a more appropriate description 
of the nonmagnetic state than what is given by the Stoner 
model. Within the collinear magnetic model, where all 
local axes of spin-quantization point in the same direc- 
tion, DLM can be treated as a binary alloy problem and 
thus described using the coherent potential approxima- 
tion (CPA)^S^ We have carried out such DLM calcu- 
lations, assuming the Cr-sublattice to be randomly oc- 
cupied by an equal number of Cr atoms with oppositely 
directed magnetic moments. The object for performing 
the DLM calculations is two-fold. If the total energy in 
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a DLM calculation is lower than the corresponding FM 
calculation, we can safely assume that the ground state 
(for the given lattice parameter and structure) is not FM, 
albeit of unknown magnetic structure. The exchange in- 
teractions in the DLM state can also be used to compute 
estimates of T c , and such estimates of T c may be consid- 
ered as estimates from above the magnetic-nonmagnetic 
transition. T c computed from exchange interactions in 
the FM state are estimates from below the transition. Of 
course, if the ground state is known to be ferromagnetic, 
then estimates of T c based on exchange interactions in 
the FM reference state are the appropriate ones to con- 
sider. 

In some cases where the FM results point to the possi- 
bility of the ground state magnetic structure being AFM 
or of complex nature, we have carried out a limited num- 
ber of AFM calculations to provide some insight into this 
problem (see section ITV Cp . 

We have computed the spin-resolved densities of states 
(DOS) for all the alloys for varying lattice parameters, 
and for both the FM and DLM configurations. The FM 
calculations show half-metallic character, due to the for- 
mation of bonding and antibonding states involving the 
t2 g orbitals of the Cr-atoms and the sp orbitals of the 
neighboring pnictogen (As, Sb) or chalcogen (S, Se, Te). 
The hybridization gap is different and takes place in dif- 
ferent energy regions in the two spin channels. The crit- 
ical values of the lattice parameters above which the FM 
calculations show half-metallic character agree well with 
those reported by Galanakis and MavropoulosJ^ The 
DOS for the alloys of the type CrX (X=Sb, S, Se, Te) 
have been presented by several other author a 12 ' 15 i 16 and 
thus will not be shown here. In Figs [1] and [2] we show the 
DOS for the mixed alloys CrAs5oSb 50 and CrAs5oSe 50 , 
for lattice parameters above and below the critical values 
for the half-metallic character. According to Galanakis 
and Mavropoulosji^ half-metallicity in ZB CrAs appears 
between the lattice parameters of 5.45 and 5.65 A. The 
latter corresponds to the lattice parameter of the GaAs 
substrate. For CrSb half-metallicity appears at a lat- 
tice parameter between 5.65 and 5.87 A. The mixed alloy 
CrAs5oSb 50 , as shown in FigJTJ is not quite half-metallic 
at the lattice parameter of 5.65 A, and fully half-metallic 
at the lattice parameter of 5.76 A. Replacing Sb with Se 
in the above alloy, i.e. for CrAssoSeso, brings the critical 
lattice parameter down slightly. As shown in Fig 121 at a 
lattice parameter of 5.65 A, CrAssoSeso is half-metallic, 
although barely so. In our calculation CrS and CrSe are 
half-metallic at a lattice parameter of 5.65 A, and not so 
at a lattice parameter of 5.55 A. CrTe is not half-metallic 
at a lattice parameter of 5.76 A, but at a lattice param- 
eter of 5.87 A. For both CrS and CrSe the critical value 
should be close to 5.65 A, and for CrTe it should be close 
to 5.87 A. 

Note that in general the half-metallic gap is larger in 
the chalcogenides than in the pnictides. This is due to 
larger Cr- moment (see section [TTTJ) for the chalcogenides, 
which results in larger exchange splitting. This explains 
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FIG. 1: Spin-resolved densities of states in ZB CrAs 5 oSb 5 o 
for lattice parameters (a) 5.55 A, (b) 5.65 A, (c) 5.76 A and 
(d) 5.87 A, respectively. 



the difference in the half-metallic gaps in Figs. [Hand [2] 
for similar lattice parameters. 

Fig. compares the total DOS of CrAs for the FM 
and DLM calculations for the equilibrium lattice param- 
eter 5.65 A. Higher DOS at the Fermi level for the DLM 
calculation, compared with the FM calculation, is an in- 
dication that the band energy is lower in the FM state. 
Indeed, as indicated in Table HI compared with the DLM 
state the total energy for ZB CrAs is lower in the FM 
state for the lattice parameters from 5.44 to 5.98 A. In 
fact, this holds for lattice parameters up to 6.62 A, show- 
ing the robustness of fcrromagnetism in CrAs over a wide 
range of the lattice parameter. This is also true for CrSb. 

In Table H] we show the variation of total energies per 
atom in Ry with the lattice parameter for CrX (X=As, 
Sb, S, Se and Te) in the DLM and FM states. The pur- 
pose of tabulating these energies is not to determine the 
bulk equilibrium lattice parameters in the ZB structure, 
as this has already been done by several authors i 12 ' 15 i 16 
Our results for equilibrium ZB phase lattice parameters 
agree with those found by Galanakis and Mavropoulos. 12 
The important point is that for CrS and CrSe at low 
values of lattice parameters the DLM energies are lower 
than the FM energies, showing clearly that the FM con- 
figuration is unstable. The result for CrS is in line with 
the observation by Zhao and Zunger^ 2 - who find ZB CrS 
to be antiferromagnetic with an equilibrium lattice pa- 
rameter of 5.37 A. As shown later ( section ITV]) , the ex- 
change coupling constants for the Cr atoms in the FM 
calculations are negative, indicating the instability of the 
ferromagnetic spin alignment. The tendency to antifcr- 
romagnetism in CrSe at compressed lattice parameters 
is also revealed in a study by Sasaioglu et alM- For 
CrTe at lower lattice parameters the FM energy is lower 
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FIG. 2: Spin-resolved densities of states in ZB CrAs 5 oSe 5 o for 
lattice parameters (a) 5.55 A, (b) 5.65 A, (c) 5.76 A and (d) 
5.87 A, respectively. 



FIG. 4: Total densities of states in CrS in the DLM and FM 
states for lattice parameters (a) 5.44 A, (b) 5.55 A, (c) 5.65 
A, and (d) 5.76 A, respectively. 
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for CrS at the lowest lattice parameter of 5.44 A the 
DOS at the Fermi level is lower in the DLM state than in 
the FM state. The deviation from ferromagnetism at low 
values of the lattice parameter for CrS, CrSe and CrTe 
is also revealed by our study of the lattice Fourier trans- 
form of the exchange interaction between the Cr atoms 
in the FM state (section IIV|) . The search for an anti- 
ferromagnetic state with lower energy is possible within 
our collinear magnetic model by enlarging the unit cell 
in various ways. We have pursued this issue to a limited 
extent, by considering 001, 111 AFM configurations for 
CrS, CrSe and CrTe at low values of the lattice parameter 
(see discussion in sectiorJlVj) . A satisfactory resolution of 
such issues is possible only by going beyond the collinear 
model. 



FIG. 3: A comparison of the total densities of states in ZB 
CrAs with lattice parameter 5.65 A for the DLM and FM 
states. 

than the DLM energy, but the exchange constants for 
the Cr-atoms in the FM calculations are still negative 
(see discussion in section IIV|) . signaling the possibility 
of the ground states in CrTe at low values of the lattice 
parameter being neither DLM nor FM. Note that in our 
discussion ground state implies the lowest energy state in 
ZB structure. For CrS, CrSe and CrTe the ground states 
at low lattice parameters can be of an antiferromagnetic 
(AFM) or complex magnetic structure. A lower total en- 
ergy may also mean a lower band energy, and in some 
cases, the latter may be reflected in a lower density of 
states at the Fermi level. This is shown in FigJU where 



III. MAGNETIC MOMENTS 

Our spin-polarized calculations for the FM reference 
states lead to local moments not only on the Cr atoms, 
but also on the other atoms (As, Sb, S, Se, and Te) as 
well as the empty spheres. Sandratskii et al£^ have dis- 
cussed the problem associated with such 'induced mo- 
ments' in case of the Heusler alloy NiMnSb and the 
hexagonal phase of MnAs. Usually such systems can be 
divided into sublattices with robust magnetic moments 
and sublattices where moment is induced under the influ- 
ence of the former. These authors argue that the treat- 
ment of the induced moments as independent variables in 
a Heisenberg Hamiltonian may lead to artificial features 
in the spin-wave spectra, but these artificial features do 
not drastically affect the calculated Curie temperatures 
of the two alloys, NiMnSb and hexagonal MnAs. Clearly, 



5 



TABLE I: Comparison of total energies in the FM and DLM states as a function of the lattice parameter for CrAs, CrSb, CrS, 
CrSe, and CrTe. Results for six lattice parameter values are shown. Calculations include six additional lattice parameters 
beyond 5.98A, reaching a maximum of 6.62 A. For all these additional lattice parameters FM energy is always lower than the 
corresponding DLM energy, indicating that ferromagnetism is favored at higher lattice parameters. 

Lattice parameter (A) 5.44 5.55 5.65 5.76 5.87 5.98 

CrAs 

DLM energy -1653.4039 -1653.4021 -1653.3995 -1653.3960 -1653.3920 -1653.3876 

FM energy -1653.4068 -1653.4055 -1653.4034 -1653.40026 -1653.3966 -1653.3922 

CrSb 

DLM energy -3762.4738 -3762.4781 -3762.4807 -3762.4821 -3762.4822 -3762.4814 

FM energy -3762.4770 -3762.4813 -3762.4841 -3762.4857 -3762.4862 -3762.4856 



CrS 

DLM energy -723.4504 -723.4468 -723.4426 -723.4381 -723.4332 -723.4280 

FM energy -723.4491 -723.4464 -723.4434 -723.4395 -723.4351 -723.4302 
CrSe 

DLM energy -1737.9146 -1737.9139 -1737.9123 -1737.9100 -1737.9071 -1737.9036 

FM energy -1737.9139 -1737.9132 -1737.9124 -1737.9110 -1737.9087 -1737.9057 
CrTe 

DLM energy -3918.9074 -3918.9124 -3918.9158 -3918.9179 -3918.9189 -3918.9188 

FM energy -3918.9078 -3918.9128 -3918.9161 -3918.9181 -3918.9194 -3918.9201 



in our case the sublattice with the robust magnetic mo- 
ment is the Cr-sublattice. Among the three other sublat- 
tices, the magnitudes of the induced moments decrease in 
the following order for the two robust ferromagnets CrAs 
and CrSb: X-sublattice (X=As, Sb), sublattice ES-1 (the 
sublattice of empty spheres that is at the same distance 
with respect to the Cr-sublattice as the X-sublattice ), 
sublattice ES-2 (sublattice of empty spheres further away 
from the Cr-sublattice). This trend is particularly valid 
for low values of the lattice parameter. The induced mo- 
ments originate from the tails of the orbitals (primarily 
d) on the nearby Cr-atoms. This is particularly true for 
the moments induced on the empty spheres. The magni- 
tudes of the induced moments on the two empty sphere 
sublattices decrease as the lattice parameter increases, 
and so do the differences in their magnitudes. The signs 
of the moments on ES-1 and ES-2 may be the same for 
small lattice parameters, but are opposite for large lattice 
parameters. The sign of the moment on the X-sublattice 
is opposite to that on the Cr-sublattice and the mag- 
nitudes of the moments on the two sublattices increase 
with increasing lattice parameters, due to decreased hy- 
bridization between Cr-rf and X-sp orbitals. Above a 
critical value of the lattice parameter, the moment per 
formula unit (f.u.) saturates at a value of 3.0 /is, as the 
half-metallic state is achieved, while the local moments 
on the Cr- and X-sublattices increase in magnitude, re- 
maining opposite in sign. The maximum ratio between 
the induced moment on X (X=As,Sb) and the moment 
on Cr is 0.18 for CrAs and 0.15 for CrSb, occurring at the 
highest lattice parameter of 6.62 A studied. The max- 



imum ratio between the induced moment on ES-1 and 
that on Cr is 0.06, occurring at the lowest lattice param- 
eter of 5.44 A studied. Magnetic moments of CrAs and 
CrSb per formula unit (f.u.) as well as the local moment 
at the Cr site are shown in Fig. [5J where we compare 
the two compounds with each other for their magnetic 
moments in the FM and DLM states. The same results 
are presented in Fig. [6j comparing the moment per f.u. in 
the FM state with the Cr local moment in the FM and 
DLM states separately for each compound. It is to be 
noted that there are no induced moments for the DLM 
reference states, i.e. the moments on the non-CR sub- 
lattices are several orders of magnitude smaller than the 
robust moment on the Cr atoms. The total moment per 
formula unit in the DLM state is zero by construction. 
The local moment on the Cr atom for the DLM reference 
state is usually less than the corresponding FM value for 
smaller lattice parameters, and larger for larger lattice 
parameters (Fig|6]). 

Similar trends in the variation of the local moment on 
Cr and the induced moments on the other sublattices for 
the FM reference states as a function of lattice parameter 
are revealed for CrX (X=S, Se, Te), except that the mo- 
ments on ES-1 are always an order of magnitude larger 
than those on ES-2. I addition, the induced moments 
on ES-2 are ~ 6 — 10 times larger than those on X sub- 
lattice for smaller values of the lattice parameter, with 
the two becoming comparable in magnitude for larger 
lattice parameters. The induced moments on ES-1 and 
X-sublattices are never larger than ~ 5% of the moment 
on the Cr atoms. The induced moments for the DLM 
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FIG. 5: Magnetic moment in ZB CrAS and CrSb as a function 
of lattice parameter in the FM and DLM states. In all cases 
studied (Figs. 0[8]), FM calculations produce 'induced mo- 
ments' on non-Cr spheres representing the X-atoms (X=As, 
Sb, S, Se, Te), and one set of empty spheres. The DLM calcu- 
lations produce no such 'induced moments', i.e., the moments 
reside on the Cr-atoms only. See text for discussion. 



reference states are several orders of magnitude smaller 
than the Cr-moments, and can be safely assumed to be 
zero. Results for ZB CrS, CrSe and CrTe are presented 
in FigsJT] and [5] The moment per f.u. reaches the sat- 
uration values of AfiB for CrS, CrSe, and CrTe in the 
half-metallic state, as discussed in detail by Galanakis 
and Mavropoulos.— The saturation values of the mo- 
ments for all these alloys (CrX, X=As, Sb, S, Se, and 
Te) satisfy the so-called "rule of 8": M = {Z tot - 8) fi B , 
where Z tot is the total number of valence electrons in the 
unit cell. The number 8 accounts for the fact that in the 
half-metallic state the bonding p bands are full, accom- 
modating 6 electrons and so is the low-lying band formed 
of the s electrons from the sp atom, accommodating 2 
electrons. The magnetic moment then comes from the 
remaining electrons filling the d states, first the e g states 
and then the ti g . The saturation value of 3/i_e/f.u., or 
the half-metallic state, appears for a larger critical lat- 
tice constant in CrSb than in CrAs. Similarly, the critical 
lattice constants for the saturation magnetic moment of 
4/is/f-u. are in increasing order for CrS, CrSe and CrTe. 
The local moment on the Cr atom can be less/more than 
the saturation value, depending on the moment induced 
on the non-Cr atoms and empty spheres. 

Fig[9]shows the variation of the magnetic moment with 
the lattice parameter for the random alloys CrAs5oX5o 
(X=Sb, S, Se, Te), where 50% of the As-sublattice is 
randomly occupied by X-atoms. The saturation moment 
per f.u. for CrAs 5 oSb 5 o in the half- metallic state is 
with the results falling between those for CrAs and CrSb 
shown in Fig|H] For CrAs5oX 50 (X=S, Se, Te), the satu- 



FIG. 6: A comparison of the moment/f.u. with the local Cr 
moment in the FM state as well as the DLM state for ZB 
CrAs (upper panel) and CrSb (lower panel). 




Lattice parameter (A) 



FIG. 7: Magnetic moments in ZB CrS, CrSe and CrTe as 
a function of lattice parameter in the FM and DLM ground 
states. 

ration moment per f.u. is 3.5/is. The local Cr- moment 
deviates from the saturation value in the half-metallic 
state, being higher than the saturation value for all lat- 
tice parameters above 6.1 A. 

From Figsj5][9] it is clear that the magnetic moment 
per formula unit is closer to the magnetic moment of the 
Cr atoms in the FM calculations than in the DLM cal- 
culations. Local Cr-moments in the DLM calculations 
are suppressed w.r.t. the FM results for low lattice pa- 
rameters and enhanced for larger lattice parameters. As 
shown in TABLE [Q the total energy of the FM state is 
lower than that of the corresponding DLM state in almost 
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FIG. 8: A comparison of the moment/f.u. with the local 
Cr moment in the FM state as well as the DLM state for 
ZB CrS(top panel), CrSe (middle panel) and CrTe (bottom 
panel) . 



all cases, except for some compressed lattice parameters 
for CrS and CrSe. However, the consideration of the 
DLM state does provide an advantage in that there are 
no associated induced moments, i.e., the DLM calcula- 
tions produce moments that reside on the robust mag- 
netic sublattice only. Mapping of the total energy on to 
a Hcisenberg Hamiltonian, therefore, does not result in 
exchange interactions involving atoms/spheres with 'in- 
duced moments' and all associated artificial/non-physical 
features referred to by Sandratskii et 
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IV. EXCHANGE INTERACTION AND CURIE 
TEMPERATURE 



Mapping onto a Heisenberg Hamiltonian and 
related issues 



Currently, most first-principles studies of the thermo- 
dynamic properties of itinerant magnetic systems pro- 
ceed via mappin g 35 ' 36 the system energy onto a classical 
Hcisenberg model: 



(i) 



where i,j are site indices, is the unit vector pointing 
along the direction of the local magnetic moment at site 
i, and J^- is the exchange interaction between the mo- 
ments at sites i and j. The validity of this procedure is 
justified on the basis of the adiabatic hypothesis- the as- 
sumption that the magnetic moment directions are 'slow 
variables' on all the characteristic electronic time scales 
relevant to the problem, and thus can be treated as clas- 



sical parameters. The energy of the system for a given 
set of magnetic moment directions is usually calculated 
via methods based on density functional theory (DFT). 

One of the most widely used mapping procedures is due 
to Liechtenstein et alM-~— It involves writing the change 
in the energy due to the deviation of a single spin from 
a reference state in an analytic form using the multiple 
scattering formalism and by appealing to the magnetic 
variant of the Andersen force theorem] 42 ' 43 The force the- 
orem, derived originally for the change of total energy 
due to a deformation in a solid, dictates that the differ- 
ences in the energies of various magnetic configurations 
can be approximated by the differences in the band en- 
ergies alonet 37 ' 38 ' 44 The energy of a magnetic excitation 
related to the rotation of a local spin-quantization direc- 
tion can be calculated from the spinor rotation of the 
ground state potential. No self-consistent calculation for 
the excited state is necessary. A second approach is based 
on the total energy calculations for a set of collinear mag- 
netic structures, and extracting the exchange parameters 
by mapping the total energies to those coming from the 
Heisenberg model given by Eq.(JTJ). Such calculations can 
be done using any of the standard DFT methods. How- 
ever, unlike the magnetic force theorem method, where 
the exchange interactions can be calculated directly for 
a given structure and between any two sites, several hy- 
pothetical magnetic configurations and sometimes large 
supercells need to be considered to obtain the values of 
a modest number of exchange interactions. In addition, 
some aspects of environment-dependence of exchange in- 
teractions are often simply ignored. The difference be- 
tween these two approaches is, in essence, the same 
as that between the generalized perturbation method 
(GPM^ii 6 . and the Connolly- Williams methodic in 
determining the effective pair interactions in ordered and 
disordered alloys. A third approach is a variant of the sec- 
ond approach, where the energies of the system in vari- 
ous magnetic configurations corresponding to spin- waves 
of different wave- vectors are calculated by employing the 
generalized Bloch theorem for spin-spirals.— The inter- 
atomic exchange interactions can be calculated by equat- 
ing these energies to the Fourier transforms of the clas- 
sical Heisenberg-model energies. This approach, known 
as the 'frozen magnon approach', is similar to the 'frozen 
phonon approach' for the study of lattice vibrations in 
solids. 

In this work, we have used the method of Liechtenstein 
et al. , which was later implemented for random magnetic 
systems by Turek et al. , using CPA and the TB-LMTO 
method^ The exchange integral in Eq. (TTJ) is given by 



Ji 



= — lim Im 

4-7T e^0+ 



tr L 



A l (z)4(z)A J (z)4 dz, (2) 



where z = E + it represents the complex energy vari- 
able, L — (/,m), and A^(z) = Pj{z) — Pf{z) } represent- 
ing the difference in the potential functions for the up 
and down spin electrons at site Y. In the present work 
g£j(z)(a =t, 4-) represents the matrix elements of the 
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FIG. 9: A comparison of the moment/f.u. with the local Cr moment in the FM state as well as the DLM state for ZB (a) 
CrAs 5 oSb 5 o, (b) CrAs 50S50, (c) CrAssoSeso, and (d)CrAs5oTe5o. Circle: moment/f.u.- FM state, diamond: local Cr-moment- 
FM state, triangle: local Cr-moment- DLM state. 



Green's function of the medium for the up and down spin 
electrons. For sub-lattices with disorder, this is a configu- 
rationally averaged Green's function, obtained via using 
the prescription of CPA. The integral in this work is per- 
formed in the complex energy plane, where the contour 
includes the Fermi energy Ep. The quantity Jij given 
by Eq. @ includes direct-, indirect-, double-exchange 
and superexchange interactions, which are often treated 
separately in model calculations. The negative sign in 
Eq.([T]) implies that positive and negative values of J,j 
are to be interpreted as representing ferromagnetic and 
antiferromagnetic interactions, respectively. 

A problem with the mapping of the total energy to a 
classical Heisenberg Hamiltonian following the approach 
of Liechtenstein et alM-~— is that it generates exchange 
interactions between sites, where one or both may carry 
induced moment (s). Of course this problematic scenario 
appears only for the FM reference states, as the DLM 
reference states do not generate induced moments. In 
the present work the Liechtenstein mapping procedure, 
applied to FM reference sates, generates exchange in- 
teractions between the Cr atoms, between Cr and other 
atoms X (X=As, Sb, S, Se, Te), and also between Cr 



atoms and the nearest empty spheres ES-1. Depending 
on the lattice parameter, this latter interaction is either 
stronger than or at least comparable to that for the Cr- 
X pairs. The exchange interactions between Cr atoms 
and the furthest empty spheres ES-2 are always about 
one or two orders of magnitude smaller than the Cr-ESl 
interactions and can be neglected. In CrAs, the ratio of 
the nearest neighbor Cr-ESl to Cr-Cr interaction varies 
from 0.2-0.25 at low lattice parameters to 0.06-0.07 at 
high values of the lattice parameter. In CrSb, these ratios 
are smaller, varying between 0.14 and 0.05. The Cr-ESl 
exchange interactions are also relatively strong in magni- 
tude in CrS, CrSe and CrTe. One important point is that 
while these interactions are positive for nearest neighbors 
for all lattice parameters, Cr-Cr nearest neighbor inter- 
action is negative for low values of lattice parameters 
in CrS and CrSe. In CrTe, this interaction changes sign 
from positive to negative and then back, as the lattice pa- 
rameter is varied in the range 5.44-6.62 A. As mentioned 
earlier, the calculation for the DLM reference states do 
not produce induced moments, and thus no exchange in- 
teractions other than those between the Cr atoms. 

Sandratskii 35 et at have discussed the case when, 
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in addition to the interaction between the strong mo- 
ments, there is one secondary, but much weaker, inter- 
action between the strong and one induced moment. In 
this case, the Curie temperature, calculated under the 
mean-field approximation (MFA), seems to be enhanced 
due to this secondary interaction, irrespective of the sign 
of the secondary interaction. In other words, the Curie 
temperature would be somewhat higher than that cal- 
culated by considering only the interaction between the 
strong moments. The corresponding results under the 
random phase approximation (RPA) have to be obtained 
by solving two equations simultaneously. One can as- 
sume that the RPA results for the Curie temperature 
follow the trends represented by the MFA results, being 
only somewhat smaller, as observed in the absence of in- 
duced moments. In our case, since there are at least two 
secondary interactions (Cr-X and Cr-ESl) to consider in 
addition to the main Cr-Cr interaction, the influence of 
these secondary interactions is definitely more complex. 

In view of the above-described situation involving sec- 
ondary interactions between Cr- and the induced mo- 
ments for the FM reference states, we have adopted the 
following strategy. Since no induced moments appear 
in calculations for the DLM reference states, the Curie 
temperature T c for these can be calculated as usual from 
the exchange interaction between the Cr-atoms, i.e. the 
strong moments. For these cases the calculation of T c 
can proceed in a straightforward manner by making use 
of the mean-field approximation (MFA) or the more ac- 
curate random- phase approximation (RPA)i^i One can 
obtain the MFA estimate of the Curie temperature from 

k B Tr A =lYi & Ci , (3) 

where the sum extends over all the neighboring shells. An 
improved description of finite-temperature magnetism is 
provided by the RPA, with T c given by 

(feD- 1 = ~ £ [ jCr ' Cr (°) - ^(q)]" 1 • (4) 

q 

Here N denotes the order of the translation^ group ap- 
plied and J Cr ' Cr (q) is the lattice Fourier transform of the 
real-space exchange integrals J^' r . It can be shown 
that T C RPA is always smaller than T c MFA .36 It h as been 
shown that the RPA Curie temperatures are usually close 
to those obtained from Monte-Carlo simulations^ As 
shown by Sandratskii et aL— the calculation of T c us- 
ing RPA is considerably more involved even for the case 
where only one secondary interaction needs to be con- 
sidered, in addition to the principal interaction between 
the strong moments. The complexity of the problem in- 
creases even for MFA, if more than one secondary inter- 
action is to be considered. The same comment applies 
to stability analysis using the lattice Fourier transform 
of the exchange interactions. The deviation of the na- 
ture of the ground state from a collinear and parallel 



alignment of the Cr moments in the FM reference states 
could be studied by examining the lattice Fourier trans- 
form of the exchange interaction between the Cr atoms: 

•^(q) = J2q ^0R Cr ex P ' ^ an tne secon dary inter- 
actions could be ignored. This is definitely not possible 
for many of our FM results, where several pairs of interac- 
tion need to be considered, and J(q) is a matrix bearing 
a complicated relationship to the energy as a function 
of the wave-vector q. Thus, in the following the results 
for T c will be presented mostly for the DLM reference 
states. For comparison, in a small number of cases we 
will present T c calculated for the FM reference states us- 
ing only the Cr-Cr exchange interactions as the input. 
Of course, this will be done with caution only for cases 
where we have reason to believe that the results are at 
least qualitatively correct. Some FM results will also be 
included towards the stability analysis based on J(q) de- 
rived from Cr-Cr interactions only. Again, this will be 
done with caution, only if the corresponding results can 
be shown to be meaningful via additional calculations . 



B. Exchange interactions for the FM and DLM 
reference states 

The Cr-Cr exchange interactions for all the alloys stud- 
ied and for both FM and DLM reference states become 
negligible as the inter-atomic distance reaches about 
three lattice parameters or, equivalently, thirty neighbor 
shells. The same applies to the Cr-X and Cr-ES inter- 
actions for the FM cases, these interactions in general 
being somewhat smaller. The Cr-Cr interactions for the 
DLM reference states are more damped compared with 
the corresponding FM results, showing less fluctuations 
in both sign and magnitude. The distance dependence 
of the exchange interactions between the Cr atoms in 
CrAs is shown in Fig[TU] for several lattice parameters. 
Although the nearest neighbor interaction is always pos- 
itive (i.e, of ferromagnetic nature), the interactions with 
more distant neighbors are sometimes antiferromagnetic. 
Such antiferromagnetic interactions are more common in 
CrAs for lower lattice parameters. With increasing lat- 
tice parameter, interactions become predominantly ferro- 
magnetic, and by the time the equilibrium lattice param- 
eter of 5.52 A is reached, antiferromagnetic interactions 
mostly disappear. We have calculated such interactions 
up to the 405th neighbor shell, which amounts to a dis- 
tance of roughly 8 lattice parameters. Although the in- 
teractions themselves are negligible around and after the 
30 th neighbor shell, their influence on the lattice sums 
continues up to about 100 neighbor shells. By about 
the neighbor 110" 1 shell (a distance of ~ 5 lattice pa- 
rameters) the interactions fall to values small enough so 
as not to have any significant effect on the calculated 
lattice Fourier transform of the exchange interaction and 
the Curie temperature (see below). It is clear from FigflUl 
that ferromagnetism in CrAs is robust and exists over a 
wide range of lattice parameters. The distance depen- 
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dence of the Cr-Cr exchange interactions in CrSb is very 
similar to that in CrAs for both FM and DLM reference 
states. 

For CrS, CrSe, and CrTe the situation is somewhat 
different. For CrS and CrSe, the FM reference states 
for some low lattice parameters yield Cr-Cr interactions 
that are antiferromagnetic even at the nearest neighbor 
separation. For CrTe, at the lowest lattice parameter 
studied (5.44 A) the nearest neighbor interaction for the 
FM reference state is ferromagnetic, but becomes antifer- 
romagnetic with increasing lattice parameter, changing 
back to ferromagnetic at higher lattice parameters. For 
all three compounds, the interactions are predominantly 
ferromagnetic at higher lattice parameters. Figs [11] and 
[T2l show the distance dependence of the exchange inter- 
actions calculated for the FM reference states in CrSe 
and CrTe, respectively, for several lattice parameters. 
Predominant nearest neighbor antiferromagnetic interac- 
tions between the Cr atoms result in negative values of 
the Curie temperature, when calculated via Eqs. © or 
((3]). These results for the Curie temperature for the FM 
reference states can be discarded as being unphysical on 
two grounds: because of the neglect of the interactions in- 
volving the induced moments and also because they point 
to the possibility that the ground state is most proba- 
bly antiferromagnetic or of complex magnetic structure. 
The antiferromagnetic Cr-Cr interactions mostly disap- 
pear, when calculated for the DLM reference states. This 
could be interpreted as being an indication that the ac- 
tual magnetic structure of the ground states for these low 
lattice parameters in case of CrS, CrSe and CrTe is closer 
to a DLM state than to an FM state. In FigQ~3]we show 
the Cr-Cr exchange interactions for the DLM reference 
states in case of CrSe for the same lattice parameters as 
those considered for FigfTTI A comparison of the two fig- 
ures shows that all interactions have moved towards be- 
coming more ferromagnetic for the DLM reference states, 
the nearest neighbor interaction for the lowest lattice pa- 
rameter staying marginally antiferromagnetic. 



C. Stability analysis via the Lattice Fourier 
transform of Cr-Cr exchange interactions 

The deviation of the nature of the ground state 
from the reference state can be studied by examining 
the lattice Fourier transform of the corresponding ex- 
change interactions between the Cr atoms: J(q) = 
^2 J^' Cr exp (iq • R). As pointed out earlier, for the 
FM reference states this procedure suffers from the draw- 
back of neglecting the effects of all other interactions in- 
volving the induced moments. For the DLM reference 
states there are no induced moments, so the relationship 
between the energy and J(q) is simpler, but a physical 
picture of the spin arrangement corresponding to a par- 
ticular wave- vector q is harder to visualize. For the FM 
reference states, if there were no moments other than 
those on the Cr atoms, a maximum in J(q) at q = 




Distance r.. (a) 



FIG. 10: (Color online) Distance dependence of the exchange 
interaction Jij between the Cr atoms in CrAs for various lat- 
tice parameters a, calculated for the FM and DLM reference 
states. The distance between the Cr atoms is given in units 
of the lattice parameter a (the same applies to Figs llU13)l . 
The main plot in Fig llOl shows the distance dependence up to 
2.25a, while the inset shows the values between 2.25a and 5a. 
Although the individual values of Jij are small beyond about 
2.25a, their cumulative effects on the total exchange constant 
and the Curie temerature cannot be neglected (see text for 
details). Comparison of the insets for the FM and DLM cases 
shows that the interactions are more damped for the DLM 
case, being at least an order of magnitude smaller for dis- 
tances beyond ~2. 25-2. 5a or 15-20 neighbor shells. Similar 
comments apply to the interactions presented in Figs llU13l 
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FIG. 11: (Color online) Distance dependence of the exchange 
interactions between the Cr atoms in CrSe for the FM refer- 
ence state. See caption of Fig llOl 
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FIG. 12: (Color online) Distance dependence of the exchange 
interactions between the Cr atoms in CrTe for the FM refer- 
ence state. See caption of FigfTUl 



FIG. 14: (Color online) Lattice Fourier transform of the Cr- 
Cr exchange interactions in CrAs, calculated for the FM and 
DLM reference states. 



\J\^CrSe: DLM 


1 1 

«-° a=5.44 A 
h-h a=5.65 A 
a=5.76 A 
«-» a=5.87 A 






\ CrSe: DLM 


a=5.98 A 
Q-a a=6.18 A 
a=6.62 A 




, , ,!,,,!,,, 



0.5 1 1.5 2 2.5 

Distance v.. (a) 



FIG. 13: (Color online) Distance dependence of the exchange 
interactions between the Cr atoms in CrSe for the DLM ref- 
erence state. See caption of Fig llOl 



would imply that the ground state is ferromagnetic with 
collinear and parallel Cr magnetic moments in all the 
unit cells. A maximum at symmetry points other than 
the T-point would imply the ground state being antifer- 
romagnetic or a spin-spiral state. A maximum at a wave- 
vector q that is not a symmetry point of the BZ would 
imply the ground state being an incommensurate spin 
spiral. The presence of induced moments and the con- 
sequent interactions involving non-Cr atoms and empty 
spheres spoil such interpretations based on J(q) derived 
from Cr-Cr interactions alone. However, the tendencies 
they reveal might still be useful. It is for this reason that 
we study the Fourier transform J(q), defined above, for 
both FM and DLM reference states. 
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FIG. 15: (Color online) Lattice Fourier transform of the Cr- 
Cr exchange interactions in CrS, calculated for the FM and 
DLM reference states. 

In FigHHwe have plotted this quantity for CrAs. The 
results for CrSb are quite similar. The maximum in J(q) 
at the T-point for all lattice parameters and for both FM 
and DLM reference states can be taken as an indication 
that the ground state magnetic structure is ferromag- 
netic for CrAs for all the lattice parameters studied. The 
same comment applies to CrSb. The apparent lack of 
smoothness in J(q) shown for the FM reference states is 
a consequence of the fact that there are other additional 
bands (involving induced moments), which are supposed 
to cross the band shown, but have not been computed. 

For CrS, CrSe, and CrTe (see Figs fTolTB) . the devia- 
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rameters the FM state has the lowest energy. For CrSe, 
AFM[111] state has the lowest energy up to the lattice 
parameter 5.66 A, as is also supported by the maximum 
of J(q) at L-point. The J(q) curves for CrSe are similar 
to those of CrS, and have not been shown. For CrTe, at 
the lowest lattice parameter of 5.44 A the lowest energy 
state is FM, as is also indicated by the maximum of </(q) 
at the T-point. For higher lattice parameters 5.65 and 
5.76A, even though the J(q) curves point to the possi- 
bility of an AFM[111] ground state, the FM state energy 
turns out to be the lowest among the configurations stud- 
ied. It could be concluded that in this case a proper rela- 
tionship between the energy and J(q), obtained without 
the neglect of the induced moments, would point to the 
ground state being FM. For these three chalcogenides, 
for lattice parameters above 5.65-5.7 A the ground state 
should be FM. 



FIG. 16: (Color online) Lattice Fourier transform of the Cr- 
Cr exchange interactions in CrTe, calculated for the FM and 
DLM reference states. 



tion of the ground state for low lattice parameters from 
the parallel arrangement of Cr moments is reflected in 
the result that the maximum moves away from the T- 
point for the FM reference states. At high values of the 
lattice parameter the maximum returns to the T-point. 
The curves for CrSe are similar to those for CrS and 
have therefore not been shown. The fact that the max- 
imum for the DLM reference sates lies at the T-point in 
most cases is again an indication that the ground state 
magnetic structure is closer to the DLM state than to 
the FM state. The conclusions based on the FM refer- 
ence state results in FigsJTS] and!6l may be suspect on 
ground of neglecting the interactions involving the in- 
duced moments. However, to explore whether they do 
carry any relevant information we have carried out ad- 
ditional calculations for the three compounds CrS, CrSe 
and CrTe for two commonly occurring antiferromagnetic 
configurations: AFM[001], AFM[111]. Note that another 
commonly occurring AFM configuration AFM[110] is not 
unique, i.e. there are several configurations that could be 
seen as an AFM [110] arrangement (see Fig3. of Ref . [53jj . 
Table 2 of Ref.[54|). The simplest among these is ac- 
tually equivalent to AFM[100]. The results for the to- 
tal energy for the two AFM calculations are shown in 
Table HJ and compared with the corresponding FM and 
DLM total energies. For CrS, the lowest energy state 
for lattice parameters 5.44 and 5.55 A is AFM [111], ex- 
actly as suggested by the maximum in J(q) appearing at 
the L-point in FigfT5]for the FM reference state and for 
these two lattice parameters. As the lattice parameter 
increases beyond 5.55 A, antiferromagnetic interactions 
diminish. For the next higher lattice parameter 5.66 A in 
Table HH the lowest energy state is DLM. This may sug- 
gest that the ground state has a complex magnetic struc- 
ture, which remains to be explored. For higher lattice pa- 



D. Curie temperatures 

We determine the Curie temperature using Eqs. 
© and flU). For the DLM reference states, 
these produce estimates of T c from above the 
ferromagnetic-s-J-paramagnetic transition, and are free 
from errors due to induced moments. However, these 
estimates are high compared with properly derived val- 
ues of T c from below the transition. The latter estimates 
would require the use of FM reference states (where the 
ground states are known to be FM) and thus a proper 
treatment of the induced moments. For CrAs and CrSb, 
the magnetic state is ferromagnetic for all the lattice pa- 
rameters considered. Hence, for the sake of comparison 
we have calculated the T c for the FM reference states 
using Eqs. ([3]) and (@| as well. According to the re- 
sults of Sandratskii et a/.— the correctly calculated T c 
values, in the presence of interactions involving all the 
induced moments, would be higher. Thus, the correct 
estimates of T c should lie somewhere between the DLM 
results and the FM results obtained with the neglect of 
the induced moments. In Fig. [T7]we show these results 
for CrAs, CrSb. We have used up to 111 shells in the eval- 
uation of Eq. ([3]) and for the lattice Fourier transform 
of J Cr ' Cr (q) in Eq. after having tested the conver- 
gence with respect to the number of shells included. The 
estimated computational error corresponding to the cho- 
sen number of shells used in these calculations is below 
±2 K. For comparison we also include the results for the 
mixed alloy CrAssoSbso, for which the calculated T c val- 
ues fall, as expected, in between those of CrAs and CrSb. 
Since RPA values are more accurate than MFA values, 
our best estimates of T c for CrAs range from somewhat 
higher than 500 K at low values of the lattice parameter, 
increasing to 1000-1100 K around the mid lattice param- 
eter range (5.75-5.9 A) and then decreasing to around 600 
K for higher lattice parameters (6.5 A and above). For 
CrSb these estimates are consistently higher than those 
for CrAs: 1100 K, 1500 K and 1200 K, respectively. The 
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TABLE II: Comparison of total energies per atom (Ry) in the FM, 
function of the lattice parameter for CrS, CrSe, and CrTe. Results for 
after the decimal, 5 only to break a tie. 



DLM, AFM[001], AFM[111], and AFM[110] states as a 
5 lattice parameter values are shown, usually to 4 places 



Lattice parameter (A) 

CrS 

DLM 

FM 

AFM[001] 

AFM[111] 

CrSe 

DLM 

FM 

AFM[001] 

AFM[111] 

CrTe 

DLM 

FM 

AFM[001] 
AFM[111] 



5.44 

-723.4504 
-723.4492 
-732.4503 
-723.4506 

-1737.91456 
-1737.9139 
-1737.9144 
-1737.91458 

-3918.9074 
-3918.9078 
-3918.9170 
-3918.9072 



5.55 

-723.4468 
-723.4465 
-723.4464 
-723.4469 

-1737.91388 
-1737.9132 
-1737.9136 
-1737.91398 

-3918.9124 
-3918.9128 
-3918.9120 
-3918.9123 



5.66 

-723.4433 
-723.4432 
-723.4420 
-723.4427 

-1737.9123 
-1737.91236 
-1737.9118 
-1737.91239 

-3918.9158 
-3918.9162 
-3918.9154 
-3918.9159 



5.76 

-723.4381 
-723.4395 
-723.4375 
-723.4382 

-1737.9100 
-1737.9110 
-1737.9094 
-1737.9101 

-3918.9159 
-3918.9181 
-3918.9174 
-3918.9180 



5.87 

-723.4332 
-723.4351 
-723.4325 
-723.4333 

-1737.9071 
-1737.9087 
-1737.9065 
-1737.9071 

-3918.9189 
-3918.9194 
-3918.9182 
-3918.9189 



estimates for CrAs are similar to those provided by Sa- 
saiolglu et alM- 

For CrS, CrSe, the results obtained with the FM refer- 
ence states would clearly be wrong, in particular, for the 
low values of the lattice parameters, for which we have 
shown the ground state to be antiferromagnetic within 
our limited search. There is a possibility that the ground 
state for certain lattice parameters might have a complex 
magnetic structure. For CrTe, even though the ground 
state appears to be ferromagnetic, there are considerable 
antiferromagnetic spin fluctuations, making the FM es- 
timates unreliable. In Fig. [18] we show the T c values for 
CrS, CrSe and CrTe for the DLM reference states. The 
values for lattice parameters for which the ground state 
has been shown to be antiferromagnetic in the preceding 
section should be discarded as being inapplicable. 

Similar results for the alloys CrAssoXso (X=S, Se and 
Te) are shown in Fig. Ql|]for DLM reference states. For 
these, the ground state is ferromagnetic for all lattice 
parameters. However, because of the neglect of the in- 
duced moments related effects, our results for the Curie 
temperatures for the FM reference states are lower than 
the properly calculated values. Thus in Fig[19]we show 
the DLM results only, which are devoid of the induced 
moment effects and provide us with estimates of T c from 
above the transition. These are expected to be somewhat 
higher than the properly computed values for FM refer- 
ence states. Thus, for these alloys the trend revealed 
in Fig. [19] for the variation of T c with lattice parame- 
ter is correct. The estimates themselves are qualitatively 
correct, albeit somewhat higher than the correct values. 
Only the RPA values are plotted in Fig. [19] which are 
more reliable than the MFA values. For comparison, we 
also show the results for the pnictides CrAs, CrSb, and 
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FIG. 17: Curie temperatures in ZB CrAs and CrSb com- 
pounds for FM and DLM reference states. For comparison, 
the results for the mixed alloy CrAssoSbso are also shown. 



CrAssoSbso, which are isoelectronic among themselves, 
but have half an electron per unit cell less than the mixed 
alloys CrAs 50 X 5 o (X=S, Se and Te) 

The differences between the results for the pnictides, 
chalcogenides and the mixed pnictide-chalcogenides can 
be summarized as follows. The pnictides, CrAs, CrSb, 
and CrAssoSbso, are strong ferromagnets at all the lat- 
tice parameters studied (5.44 A- 6.62 A). In the DLM de- 
scription, their T c stays more or less constant (apart from 
a minor increase) as the lattice parameter increases from 
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FIG. 18: Curie temperatures in ZB CrS, CrSe and CrTe calcu- 
lated for DLM reference states. For CrS and CrSe the ground 
state is antiferromagnetic (see text for details) for some low 
values of the lattice parameter. The values shown for those 
lattice parameters should be discarded. 



1250 - 



1000 - 



750 



500 



DLM 




•— • CrAi - RPA 




♦— ♦ CrSb - RPA 




A- -A CrAs SQ Sb^ 


RPA 


W- V CrAs g(> S g|( - 


RPA 


>- > CrAs ;u Se ;0 


- RPA 


0-0 CrAs 5o Te 5 o 


RPA 



5.4 5.6 5.8 6 6.2 6.4 

Lattice parameter (A) 



6.6 



FIG. 19: Variation of Curie temperature as a function of lat- 
tice parameters in ZB CrAssoXso alloys with X=S, Se and Te. 
For comparison, the results for CrAs, CrSb and CrAssoSbso 
are also shown. All results shown are for DLM reference 
states, and as such, should be considered as upper limits for 



5.4 /AA to 6.1 A, and then decreases beyond (Figs. [T7] 
and[l9j The chalcogenides are antiferromagnetic or have 
complex magnetic structure for low lattice parameters. 
In the DLM description, their T c in the ferromagnetic 
state increases and then becomes more or less constant as 
the lattice parameter increases (Fig. [L8|) . The mixed al- 
loys CrAs5oX 50 (X=S, Se, Te) are ferromagnetic at all the 
lattice parameters studied. In the DLM description, their 
T c rises and then falls as the lattice parameter is increased 



from 5.44 Ato 6.62 A. A comparison of the results pre- 
sented in Figs. [L71IT91 shows that large changes in T c take 
place by changing the number of carriers. Changes due 
to isoelectronic doping are small compared with changes 
brought about by changing carrier concentration. 



V. SUMMARY OF RESULTS 

Our ab initio studies of the electronic structure, mag- 
netic moments, exchange interactions and Curie tem- 
peratures in ZB CrX (X=As, Sb, S, Se and Te) and 
CrAs 5 oX 5 o (X=Sb, S, Se and Te) reveal that half- 
metallicity in these alloys is maintained over a wide range 
of lattice parameters. The results for the exchange inter- 
action and the Curie temperature show that these alloys 
have relatively high Curie temperatures, i.e. room tem- 
perature and above. The exceptions occur for the al- 
loys involving S, Se and Te at some low values of lattice 
parameters, where significant inter-atomic antiferromag- 
netic exchange interactions indicate ground states to be 
either antiferromagnetic or of complex magnetic nature. 
A comparison of total energies for the FM, DLM, and 
two ZB antiferromagnetic configurations (AFM[001] and 
AFMflll]) show the lowest energy configuration to be 
AFMflll] for CrS and CrSe for compressed lattice pa- 
rameters (TablfJTT|. The possibility of AFM ground states 
for compressed lattice parameters for CrS was noted by 
Zhao and Zunger 22 and for CrSe by Sasioglu et al.M Our 
search for the antiferromagnetic ground states is more 
thorough than what was reported in these two studies. 
An extensive study of several antiferromagnetic configu- 
rations as well as ferrimagnetic and more complex mag- 
netic structures for CrS, CrSe and CrTe is currently un- 
derway. 

The mixed pnictide-chalcogenide alloys CrAssoXso 
(X= S, Se, Te) do not show any tendency to antiferro- 
magnetic spin fluctuations for the entire range of the lat- 
tice parameter studied. Presumably the pnictogens sup- 
press antiferromagnetic tendencies. Such alloys may play 
an important role in fabricating stable ZB half-metallic 
materials, as the concentration of the pnictogens and 
the chalcogens may be varied to achieve lattice-matching 
with a given substrate. As long as the concentration of 
As or Sb is higher than the chalcogen concentration, half- 
metallic ferromagnetic state can be achieved. There is a 
large variation in the Curie temperature of these alloys 
(Fig. [T^|) as the lattice parameter varies from the low 
(~ 5.4 A) to the mid (~ 6.1 A) range of the lattice pa- 
rameters studied. This variation is much smaller for the 
isoelectronic alloys CrAs, CrSb and CrAssoSbso (Fig. [T7)) 
over this range of lattice parameters. Note that most II- 
VI and III-V ZB semiconductors have lattice parameters 
in this range. Large changes in T c can be brought about 
by changing the carrier concentrations. The pnictides in 
general have a higher T c than the chalcogenides. 

Our results for the Curie temperature, the lattice 
Fourier transform of the exchange interactions, and the 
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resulting stability analysis are based on the exchange in- 
teractions between the Cr atoms only. For the FM ref- 
erence states this causes some errors due to the neglect 
of the effects of the induced moments. The DLM results 
are free from such errors. It is expected that the present 



study will provide both qualitative and quantitative guid- 
ance to experimentalists in the field. 
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